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ABSTRACT 

We make a detailed analysis of the exact relativistic core-shell models recently proposed 
to describe a black hole or neutron star surrounded by an axially symmetric, hollow halo of 
matter, and in a seminal sense also galaxies since there are massive shell-like structures - as 
for example rings and shells - surrounding many of them and also evidences for many galactic 
nuclei hiding black holes. We discuss the unicity of the models in relation to their analyticity 
at the black hole horizon and to the full elimination of axial (conical) singularities. We also 
consider Newtonian and linearized core-shell models, on their own to account for dust shells 
and rings around galaxies and supernova and star remnants around their centers, and also as 
limiting cases of the corresponding relativistic models to gain physical insight. 

Secondly, these models are generic enough to numerically study the role played by the 
presence/lack of discrete reflection symmetries about planes, i. e. the presence/lack of equatorial 
planes, in the chaotic behavior of the orbits. This is to be contrasted with the almost universal 
acceptance of reflection symmetries as default assumptions in galactic modeling. We also 
compare the related effects if we change a true central black hole by a Newtonian central mass. 
Our main numerical findings are: 

1- The breakdown of the reflection symmetry about the equatorial plane in both Newtonian 
and relativistic core-shell models does i) enhance in a significant way the chaoticity of orbits 
in reflection symmetric oblate shell models and ii) inhibit significantly also the occurrence of 
chaos in reflection symmetric prolate shell models. In particular, in the prolate case the lack of 
the reflection symmetry provides the phase space with a robust family of regular orbits that is 
otherwise not found at higher energies. 

2- The relative extents of the chaotic regions in the relativistic cases (i. e. with a true central 
black hole) are significantly larger than in the corresponding Newtonian ones (which have just a 
— 1/r central potential) . 



Subject headings: black holes — galaxies: structure — galaxies: rings, shells — circumstellar 
matter: rings, shells — dark matter: hollow halos — stars: stellar dynamics 



1. Introduction 

We make a more detailed study, both analytical and numerical, of some exact relativistic solutions 
we recently proposed to describe static, axially symmetric massive core-shell systems (Vieira and Letelier 
1996a, 1997, hereafter VL1 and VL2, respectively). We add in this way to the efforts of modeling many 
situations of interest in astrophysics involving massive shell-like structures around centers, as for example 



- 2 - 



black holes or neutron stars surrounded by massive shell and ring remnants. A nice illustration of this 
possibility is offered by the famous Supernova 1987A plus its physical rings, see e. g. Panagia et al. (1996), 
Meyer (1997) and Chevalier (1997). At least in a seminal sense, the model could also describe galaxies 
since there are many of them exhibiting massive rings (Sackett and Sparke 1990, Arnaboldi et al. 1993, 
Reshetnikov and Sotnikova 1997) and shells (Malin and Carter 1983, Quinn 1984, Dupraz and Combes 
1987, Barnes and Hernquist 1992), while many others possibly have galactic nuclei hiding black holes (see 
Kormendy and Richstone 1995 for a review). 

Specifically, we implement here a monopolar core (a black hole in the relativistic case) plus an exterior 
shell of dipoles, quadrupoles and octopoles. Obviously, these multipoles are shell- like Legendre expansions 
i. e. their corresponding terms increase with the distance in the intermediate vacuum between the core and 
the shell. Beyond its own applicability, the model also points to the possibility of a realistic description of 
any axially symmetric relativistic core-shell configuration whose approximation as a static system should 
be valid at some useful time scale, via a systematic multipolar Legendre expansion. In fact, a further step 
in this program was lastly started in Letelier and Vieira (1997) by giving stationarity (rotation) to the 
relativistic case of a monopolar core plus a purely dipolar shell. This is an important improvement since 
real celestial objects do rotate. The present analysis will be extended to rotating core-shell models in a 
forthcoming contribution. 

Additionally, we consider the Newtonian counterpart of the relativistic model above, which was only 
sketched in the previous works, and show that they could describe core-shell systems "per se" interesting in 
astronomy. So, we can have massive circumstellar dust shells around certain types of stars as byproducts of 
their death, see e. g. Barlow et al. (1994) for an example of (irregular) circumstclar shells around luminous 
blue variables and a review by Groenewegen et al. (1998) about dust shells around carbon Mira variables. 
Another potential application for core-shell models is to the more speculative possibility of hollow galactic 
halos of dark matter made of neutrinos recently considered by Ralston and Smith (1991), Madsen (1991) 
and Barnes (1993). Moreover, we will see that the study of the relativistic core-shell model in parallel to its 
Newtonian counterpart is useful also to clarify the physical content of the former one. 

After Poincare (1957) and the KAM theory (after Kolmogorov 1954, Arnol'd 1963a, 1963b and Moser 
1967) it became well established that non-integrability and hence chaos is a general rather than exceptional 
manifestation in the context of dynamical systems (modeling or not physical situations) (see Berry 1977). 
Given this ubiquitous fact, an important issue in astronomical modeling is to study in which extent in phase 
space chaoticity rises in models that are relevant to describe real systems and what are its consequences. 
For example, Binney (1982a) discusses the difficulties of constructing stationary self-consistent models 
when a significant fraction of orbits are irregular since they may not obey neither Vlasov's nor Jean's 
equations. Then, it is remarkable that a wide class of fully integrable potentials, the so called Stackel 
potentials (Lynden-Bell 1962, de Zeeuw 1985 and de Zeeuw et al. 1986), are feasible starting points to 
describe, by themselves or by adding perturbations, real disc, elliptical or even triaxial galaxies. For more 
realistic tridimensional models there are evidences for rounder and smoother mass distributions generating 
only relatively small fractions of chaotic orbits (Schwarzschild 1979, Binney and Tremaine 1987 and Evans 
et al. 1997), while flattening and/or sharpening of the mass distributions, for example through increasing 
triaxiality and/or putting cusps and central masses to mimic black holes, tends to increase the chaoticity 
and force us to take it into account (Gerhard and Binney 1985, Schwarzschild 1993, Merrit and Fridman 
1996, Norman et al. 1996, Merrit 1997 and Valluri and Merrit 1998). On the other hand, the emergence of 
chaos in two-dimensional models has more loose correlations with morphological aspects, see for example 
Binney and Spergel (1982), Richstone (1982), Binney (1982b), Gerhard (1985) and Sridhar and Touma 
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(1997). 

We address in this work also a numerical study about the chaoticity of orbits trapped in the bound 
gravitational zones between the core and the external massive shell. Beyond its own applicability as 
seen above, we take advantage of models that are generic within the class of axially symmetric core-shell 
distributions and moreover offer exact relativistic and Newtonian counterparts to be compared, in order 
to achieve some understanding about chaoticity related to two aspects: firstly, to the role of reflection 
symmetries shared by almost all models in astronomy and astrophysics, and secondly to the consequences 
of treating exact central black holes as just Newtonian central masses. 

In the first part of this article (section ||) we present analytical results concerning some properties 
like unicity and analyticity of the relativistic solutions themselves, in connection with both the linearized 
core-shell solutions and the corresponding Newtonian models taken as limiting cases. 



In section 2.1 we enlarge the relativistic solution presented in VL1 for a monopolar core plus a shell of 
quadrupoles and octopoles in order to include a dipolar shell component. This complete solution was only 
outlined in VL2. In this later reference it was emphasized that the case of a core plus purely dipolar shell is 
physically non-trivial already in the Newtonian gravity, moreover, it was shown that the Newtonian dipolar 
case is integrable, in contrast with its chaotic relativistic counterpart, which allowed us to characterize 
the chaotic behavior in this case as an intrinsic general relativistic effect. Here, we take advantage of 
this striking difference of the dipolar case in order to reinforce the numerical conclusions presented in 
section || concerning the non- trivial dynamical role played by the presence/absence of the discrete reflection 
symmetry about the equatorial plane on the chaoticity of the orbits. 



In section 2.2 we proceed to a full elimination of axial (conical) singularities (hereafter CSs) from the 
relativistic model, which was only partially accomplished in VL1. In fact, CSs are not globally removable 
from neither static nor stationary many-body relativistic solutions since they are self-consistently demanded 
by Einstein's equations to strut the otherwise unstable configuration against its own gravity (see Robertson 
and Noonan 1968 and Letelier and Oliveira 1998 for a recent review). So, in the present case the most we 
can do is really move all them outside the shell, thus obtaining a true vacuum in the intermediate space 
between the core and the shell. We also discuss in this section the conditions for the Kruskal-type analyticity 
of the solutions at the horizon of the central black hole and use both aspects to find the conditions under 
which the unicity of the models is assured. 



In section 2J5 we present the linearization of the relativistic model in the multipole strengths via the so 
called Regge- Wheeler (RW) formalism. In particular, we see that our model exemplifies the fact that CSs 
survive to the linearization process, their presence in the intermediate vacuum being in fact an obstruction 
to the application of the RW formalism. This does reinforce the interpretation of CSs as a kind of singular 
matter distribution necessary to the dynamical consistency of static and stationary relativistic models. 
However being intrinsecally relativistic manifestations, an adiabatic treatment of CSs in the Newtonian 
limit allows for estimates concerning the emission rate of gravitational waves by two coalescing rotating 
black holes (Araujo et al. 1998). 



In section 2.4 we discuss Newtonian core-shell models, on their own as well as limiting cases of the 
corresponding relativistic ones. We put the Henon-Heiles-like structure of the relativistic models in the 
appropriate astronomical context and stress the physical content of the various terms present in them, in 
particular we see that the apparently naive constant relativistic solution does in fact hide a relativistic 
homoeoid (see Chandrasekhar 1987 for Newtonian homoeoids). 
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The numerical part of this article is presented in section ^ and deals with a study of the chaoticity of 
orbits trapped in the intermediate vacuum between the core and the shell. Specifically, we shall explore 
in this section the fact that, due to its generality, axially symmetric core-shell expansions are particularly 
suitable to study the dynamical role played by a discrete symmetry of the model — namely the reflection 
symmetry about an equatorial plane, on the chaoticity of orbits. The reflection symmetry is present in the 
model if and only if the shell is of even-type, i. e. iff 2 2n -poles, n — 0, 1, 2, 3, occur in the shell expansion, 
and it is broken if we add any odd, 2 2n+1 -poles to the expansion. The hypothesis of reflection symmetry 
is a widely spread assumption in astronomical modeling and its twofold justification lies firstly in that too 
many real celestial objects to be modeled, for example stars themselves, star clusters and galaxies seem 
in fact very symmetric with respect to a middle plane. Another reason for this symmetry assumption is 
that the resulting model gets strongly simplified in both analytical and numerical aspects. Of course, we 
do not expect the reflection symmetry be realized in nature exactly, so it is relevant to search for possible 
detectable dynamical effects arising from deviations of that symmetry. 

A strong additional motivation for this numerical study is to compare from the point of view of the 
chaotic behavior what happens when the presence of true central black holes are simplified by reducing 
them to Newtonian central masses. This is a commom practice in current modeling, as for example in 
Gerhard and Binney (1985), Sridhar and Touma (1997) and Valluri and Merrit (1998). 

Surprisingly enough, the numerical findings show marked effects on the chaoticity of the orbits in the 
intermediate vacuum between the core and the shell linked to the presence/lack of the reflection symmetry 
in the relativistic as well as Newtonian models. We also find strong quantitative differences in the chaoticity 
manifested in both relativistic and Newtonian cases. 

Finally, we present the conclusions with some discussion and prospects in section H. 



2. Analytical results: exact, linearized and Newtonian core-shell models 

2.1. The exact model 

We deal here with static, axially symmetric models, for which the Weyl coordinates (t,p,z,tfi) are the 
starting point: 

ds 2 = e 2v dt 2 - e 2 ^- 2v [dz 2 + dp 2 } - e- 2v p 2 d4> 2 , (1) 

where v and 7 are functions of p and z only. Except where units are explicitly required (particularly in 
section 



2.4), we use non-dimensional variables: s «-> s/L (idem for p and z), and t <-> ct/L where c is the 
light velocity in vacuum and L is some convenient unit of length. Einstein's equations in absence of matter 
(vacuum) reduce in this case to the usual Laplace equation for v 

v, PP + - f, P + v,zz = 0, (2) 

and the quadrature 

d l = p[(^, P ) 2 - (v,z) 2 } dp + 2pu jP v >z dz (3) 

for 7. 

Before proceeding, we mention that the Weyl coordinates should be considered somewhat deceiving if 
we insist in naively transferring their image contents about mass configurations to the Newtonian common 
sense. For example, the spherical shape of the horizon of a black hole is compressed into a bar in Weyl's 
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coordinates, the interior portion of the black hole space-time being wholly discarded from the portrait. 
Apart from the well known fact that the relativistic context is not cast in a straight relation with the 
Newtonian one in strong regimes, they are mathematically sound as will be clear later when we will compare 
both versions of the core-shell models. 

We pass to the prolate spheroidal coordinates (t, u, v, (j>), which have a direct link with the "spherical" 
ones (t, R, 0, 4>) (R <-» R/L) we will use later, 

U=R-l= ±[vV + (z + 1)2 + y/p2 + (z- 1)2], U >1, 

v = cos6 = ^[Vp 2 + (z + I) 2 - vV + {z- l) 2 ], -1 < v < 1, (4) 

or, in terms of p, z, 

p = v/(m 2 - i)(i - v 2) = y/R(R-2)sin6, R>2, 

z = uv = (R— l)cos6. (5) 
Eqs. (Q) and (|^) are written in terms of u, v as 

[(u 2 - l)i/„] iU + [(1 - v 2 )u v ], v = 0, (6) 



l.u = n[u(u 2 - l)(v. u ) 2 - u(l - v 2 ){v, v ) 2 - 2v(u 2 - l)v u v v \, 

"f,v = ^\[v{u 2 - l){v, u ) 2 - v(l - v 2 ){u v f + 2u(l - v 2 )u u u v \. (7) 
u — v z 

In the prolate spheroidal coordinates u, v, Laplace equation (^) can be separated and solved in terms 
of standard Legendre polynomials Qg, Pi (see for example Moon and Spencer 1988, and, for relativistic 
applications, Reina and Treves 1976): 

oo 

u(u,v) = Y^[aiQi(u) + biPi(u)}[c e Qi(v) + d e Pi(v)}. (8) 
o 

The particular solution picked out from the general one above is determined by the matter distribution 
whose model is wanted. We are interested here in monopolar core plus shell-type models, so we are guided 
by the Newtonian case (to be detailed in section |2.4| ) to the specific solution of the form (up to third order) 

2v = 2^o - 2kQ (u) + 2VP 1 (u)P 1 (v) + (4/3) QP 2 (u)P 2 (v) + (A/5)OP 3 (u)P 3 (v). (9) 

va is an integration constant and Qq describes the monopolar core, which either reduces to a black hole if 
we put k = 1 and identify 2L with the Schwarzschild radius of the core (in which case (t, R, 6, 4>) above 
are Schwarzschild's coordinates), or can be switched off by simply putting k — 0. The remaining terms 
correspond to the multipoles originated from the exterior shell of matter: dipole T>, quadrupole Q and 
octopole O (note the opposite sign convention for dipoles in this definition with respect to that appearing 
in VL2). The nontrivial character of shell dipoles in both Newtonian and general relativistic theories of 
gravity has been anticipated in VL2 while shells made of quadrupoles plus octopoles were considered in 
VL1. If non-monopolar cores should be considered, we should add to (||) terms of the form Qi(u)Pt(v), 
£> 1. 
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Much more sophisticated multipolar relativistic treatments for core, shell and core-shell models are 
considered by Thorne (1980), Zhang (1986) and Suen (1986a, 1986b), respectively, in terms of higher order 
systematic expansions of the metric in de Donder coordinates. Our approach is more modest at this stage 
and inspired only on the zeroth order, Newtonian limit of the full relativistic situation. 

The odd multipoles of v in © (V and O in the present case) break the reflection symmetry about the 



plane z = 0, and we shall show in section 2.2 that we need of both V and O being simultaneously either 
present or absent if we want to rule out conical singularities of the intermediate vacuum between the core 
and the shell. We rewrite @ as 

2i/ = Klog(s=i) + P( U ,«) 1 

P = 2v + 2Vuv + ±Q(3u 2 -l)(3v 2 - 1) + \Ouv{5u 2 - 3)(5v 2 -3). 1 ' 



After integrating (^) we obtain the corresponding solution for 



2 1 = K 2 \og(^) + Q(u,v), 

Q = 270 + 7D + 7Q + 70 + JDQ + 1DO + 7QO, 

713 = AnVv - V 2 [u 2 (l - v 2 ) + v 2 ], 

7 Q = -4kQu(1 - v 2 ) 

+ (1/2)Q> 4 (1 - v 2 ){\ - 9v 2 ) - 2u 2 (l - v 2 ){\ - 5v 2 ) - v 2 (2 - v 2 )}, 
7o = -(2/5)kOv[5(3u 2 - 1)(1 - v 2 ) - 4] 

+ (3/100)O 2 [-25u 6 (l-w 2 )(5w 2 + 2w-l)(5« 2 -2u-l) 1 ' 

+ 15u 4 (l - w 2 )(65w 4 - 40w 2 + 3) - 3u 2 (l - v 2 )(25v 2 - 3)(5w 2 - 3) 

-u 2 (25u 4 -45u 2 + 27)], 
lDQ = -4PQ™(u 2 - 1)(1 - v 2 ), 

1DO = (3/10)2?O[u 2 (5u 2 - 6)(1 - v 2 )(l - bv 2 ) + v 2 {5v 2 - 6)], 

7QO = ~{<6/b)QOuv{u 2 - 1)(1 - w 2 )[(5u 2 - l)(3v 2 - 1) + 2(1 - v 2 )]. 

This is the complete solution, which was outlined in VL2. The additional integration constant 70 will be 
fixed in the next section in connection with the elimination of conical singularities from the intermediate 
vacuum, while the integration constant vq will prove to be necessary in assuring analyticity to the full 
core-shell solution at the horizon in the absence of conical singularities. The terms proportional to k in jd, 
7q and 70 represent nonlinear interactions between the black hole and the external shell. 



2.2. Unicity and smoothness requirements 

We shall be interested in two classes of singularities of the Weyl spacetime. The first one are the 
strong singularities that are located in the points wherein the scalar polynomial invariants of the curvature 
tensor blow up. For a static axially symmetric spacetime solution to the vacuum field equations we have 
only two non-vanishing invariants. They are (Carminati and McLenaghan 1991) wi = \C a bcdC ahcd and 
u>2 = ^jgC-^C^C^f where C abcd is the Weyl trace-free tensor, which for vacuum solutions coincides with 



-7- 



the Riemann curvature tensor. After some algebraic manipulations, they reduce to 

wi = 2n 2 \ia{y 2 z + p 2 v 2 p v 2 z - v, p /p - pv, P o) + v 2 p (\ + 2fr/ iP + 3pv iP ) 

+ P 2 {v% + v %) + v, P z(6v, P v,z + v, P z - 2{v,z7, P + v,p1, z )) 

+ f,p P (37,p(l - 2 P V , P )/P - v ,zz + 1v,zl,z)\, 
w 2 = 3K 3 p- 2 (u p /p - o)[o(Zp 3 v 2 p v 2 z + p(l - 3pv p )(J + 2u p ) (12) 

- v 2 z (7 P + 3^(1 - 2pu p )) + p3{v% + v %) 

+ rv iPZ (2v iZ (3i> tP (l - pv, p ) + pv 2 z ) + v pz ) 

+ "ppiHpi 1 - 2 PKp) + 4rV p - pv zz )} 

where a = v\ + v 2 r an n = cxp 2 (y — 7) . In the present case, these scalars are singular on the position of the 
attraction center and will be also singular in some directions of the spatial infinity (the specific directions 
will depend on the signs of T>, Q, etc.), in other words they are singular in the position of the sources of the 
gravitational field. 

There is another class of singularities that do not show up in a simple way in the curvature, the so 
called conical singularities (CSs) (Sokolov and Starobinskii 1977). These singularities are in some sense 
like the distributions related to a low dimensional Newtonian potential, e.g., for an infinite massive wire 
the potential is <f> = 2A In p, we have that the Laplace equation, the analog to the curvature, gives us 
V 2 = for p 7^ while for p = we need to use some global property, for example Gauss theorem, 
to get W 2 (p — 47rA(5(p). To be more precise, let us consider a conical surface, z = ap, embedded in 
the usual Euclidean three dimensional space ds 2 = dp 2 + p 2 dip 2 + dz 2 , so that we have on the cone 
ds 2 = (1 + a 2 )dp 2 + p 2 dip 2 = dp 2 + (3 2 p 2 d(p 2 with p = vT+ a 2 p and (3 2 = 1/(1 + a 2 ). Note that the 
coordinate range of p as well as ip are the usual ones. The ratio between the arc of the circumference and 
its radius is 2[3tt in this case. If we compute the Riemann tensor for this last metric we find that R pippv = 
for p ^ 0. By using the Gauss-Bonnet theorem in this case, we can pick up the curvature singularity as 
being of the form R pvpv oc S(p). 

It is well known that CSs arise when we consider Weyl's solutions near the symmetry axis (Robertson 
and Noonan 1968 and Letelier and Olivcira 1998a). They are interpreted as a geometric consequence of the 
presence of some kind of "strut" , necessary to consistently prevent any static non-spherically symmetric 
model from collapsing due to self-gravity. In this sense, the uniqueness of the core-type, one-body solution 
for the (static) spherically symmetric Einstein's equations (namely, the Schwarzschild one) should relate 
to the fact that it is impossible to attach struts in a perfectly round metric. For the Wcyl solutions, we 
can always consider a small disk centered in the symmetry axis (t = t 0} z = zo, p = e, < (p < 2ir) and 
impose that the ratio between the circumference and its radius equals 2tt. This is the condition to eliminate 
CSs from the intermediate vacuum (then attaching them to the infinity), which amounts to impose on the 
function 7 the well known conditions of elementary flatness 

27 \p=0,\z\>l= 27 \u=\z\,v=±l= 0. (13) 

These two conditions fix the constant 70 and impose an additional constraint on the odd shell multipoles 
(D and O here) in the presence of the black hole: 

^o-V 2 - 1 -Q 2 -^- Q 2 -^VO = 0, (14) 



n[V+\o] =0. 



(15) 
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The condition ( |l4| ) fixing 70 is necessary in all cases to rule out the conical singularities. From (|T^) wc 
see that the former is also a sufficient condition in two cases: i) the shell of dust is left alone by switching 
off the core (k = 0), and ii) in the presence of the black hole (k — 1), the shell is made only of even-type 
multipoles (Q here). We also see that it is possible to eliminate conical singularities for a single or pure 
shell component only if it is of the even type. Then, it follows that the core-shell dipole solution presented 
in VL1, as well as those with ^ presented in VL2, all have conical singularities since they do satisfy 
(U) but not (|l|) . If we include the necessary condition ([l4]) in (|ll]), we have the following for 7: 

27 = « 2 l g(^fi) + Q( u ,«), 

Q = JD + 7Q + 70 + JDQ + JDO + 1QO, 

7D = AkDv - V 2 (u 2 - 1)(1 - v 2 ), 

Iq = -4kQu(1 - v 2 ) - (l/2)Q 2 (u 2 - 1)(1 - v 2 )[u 2 (9v 2 - 1) + 1 - v 2 ], 
70 = -(2/5)nOv[b(3u 2 - 1)(1 - v 2 ) - 4] 

-(3/100)O 2 (u 2 -l)(l-v 2 )[(25u i -20u 2 + l)(25v i -Uv 2 + l) 1 ' 

+ 30u 2 v 2 (5v 2 - 1) + 6(1 - v 2 )), 
7dq = -AVQuv{u 2 - 1)(1 - v 2 ), 
7do - (S/l0)VO(u 2 - l)(5u 2 - 1)(1 - v 2 )(l - bv 2 ), 
lQO = -{<o/b)QOuv{u 2 - 1)(1 - v 2 )[{bu 2 - l)(3w 2 - 1) + 2(1 - v 2 )]. 

Here, it is explicitly seen that only the first terms in 723 and 70, respectively, do not vanish identically 
at v = ±1. They cancel one another at those points only if we use the additional condition ([[5]) in the 
equation above, which amounts to put there either n — or V + (2/5)0 = plus k = 1. In the later case, 7 
lastly becomes 

2 1 = log(^) + Q(u,v), 

Q — 7Q + IDO + 1QDO 1 

7Q = -4Qu(l - v 2 ) - (l/2)Q 2 (u 2 - 1)(1 - v 2 )[u 2 (9v 2 - 1) + 1 - v 2 }, 
lDO = -20v{iu 2 - 1)(1 - v 2 ) - (l/4)0 2 (w 2 - 1)(1 - v 2 )x 

(75u 4 v 4 - A2u i v 2 - A2u 2 v A + 18u 2 v 2 + 3u 4 + 3v 4 + 1), 
Iqdo = -2QOuv(u 2 - 1)(1 - v 2 )(3u 2 - l)(3v 2 - 1), 

remembering that now D is present only through 0. 

It remains an arbitrariness in the full solution, namely the constant vq. To show that it plays a 
non-trivial role in assuring analyticity to the solution at the black hole horizon, we start by writing the 
solution in the Schwarzschild coordinates (t, R, 9, </>): 

ds 2 = (1 - -)e p dt 2 - e Q " p [(l - -)" W + R 2 d6 2 } - e~ p R 2 sin 2 6d0 2 , (18) 
R R 

with P = P(u = R — l,v = cosO) and Q — Q(u = R — 1, v = cosQ) given respectively by (|Io| ) and (|ll| ) 
with k — 1. At first sight we should expect that the singularity of this metric at the horizon remains only 
a coordinate defect since all metric functions here do differ from the corresponding Schwarzschild ones by 
well-behaved exponentials. To see what really happens there, we go to Kruskal coordinates (V,U,6,<j>), 
defined as 

{R-2)e R ' 2 = U 2 -V 2 , 

t = 21og(^), (19) 

in order to eliminate the usual horizon divergence coming from the factor (1 — 2/R)^ 1 . The new, Kruskal 
components g' (x') are obtained from the old, Schwarzschild ones g a p(x) via 

= IS J^ (X) - (20) 



(17) 
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The line element ( |18| ) reads in Kruskal's coordinates as 

ds 2 = F 2 [-(l - H ■ U 2 )dU 2 + (l + H- V 2 )dV 2 ] 

-eQ- p R 2 d0 2 -e- p R 2 sin 2 9d0 2 , { ' 

where R (and hence u = R — 1) is implicitly, analytically given in terms of U, V by the first relation in (JTgj) , 
v = cos9, and F 2 and H are defined by 

F 2 = M e -R/2 e P 

R ■ (22) 

H={l~eQ- 2P ){U 2 -V 2 )- 1 . [ ' 



An inspection of ( |21| ) and (22) clearly shows that this metric is analytic at the horizon (u = R — 1 = 1 
or U = ±V) if and only if H is analytic there. If we define u = R — 1 = 1 + e, from which follows 
that U 2 — V 2 — e e^ 1+t l 2 \ this amounts just to impose the finiteness of lim c ^o H . In fact, we see that 
Q — 2P = Co + f(cos6)e + 0(e 2 ), with Co a constant depending only upon v , 70, V, Q and O, and / being 
a function only of cos9, so that the limit 

lim H = lim{e- (1+e/2) _ e c °f(cos9) + 0(e)]} (23) 

is finite (and equals to —e~ 1 f(cos9)) if and only if the condition Co = holds identically. This additional 
constraint on z/ and 70 reads 

-2(2^o + f Q) + [2 7 o ~ V 2 - \Q 2 - ^O 2 - ±VO] = 0, (24) 

which fixes uq. We see that the analyticity at the horizon and the non-existence of conical singularities are 
independent conditions, however if the necessary condition for the absence of conical singularities ( |l4| ) also 
holds, ( pi| ) reduces to the following first order relation 

2^o + ^Q = 0. (25) 

In any case both conditions taken jointly make the solution to be unique. 

It is worth to emphasize that assuring analyticity to the metric at the horizon has nothing to do with 
removing neither strong nor conical singularities from curvature invariants, since the vacuum between 
the core and the shell are free from strong singularities while the conical ones are removable from there 
through conditions (|l4j) and (|l5|) above. Secondly, however Kruskal coordinates are historically linked to 
the question of geodesic completness, this issue is not necessary to our purposes. In fact, we use here the 
analyticity of the Kruskal coordinates near the horizon of the central black hole only to assure that the 
shell is a truly controllable perturbation of the black hole in that region if expanded at any truncation 
order in the shell parameters, in the same sense considered by Vishveshwara (1970) for linear perturbations. 
Relations (|2~i| ) and (|25| ) show that analyticity near the horizon is a nontrivial property of the relativistic 
shells, having to be forced into the solution to assure its analytic behavior there. 



2.3. The linearized model: conical singularities and the Regge- Wheeler formalism 



We consider now the linearization of the solution with respect to the shell parameters, i. e. , we 
consider the shell as a small perturbation (e. g. formed by dust) of the central black hole geometry. Due to 
the spherical symmetry of the Schwarzschild background, it is assumed that any linear perturbation could 
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be expanded in spherical harmonics, an approach due firstly to Regge and Wheeler (1957) to study the 
stability of black holes under small perturbations. We shall see that this is not true if conical singularities 
are present, in spite of the perturbation we are dealing with being a perfectly linear one. 

We firstly summarize the Regge- Wheeler approach. The metric is expanded as g^ v = g^ lv + eh )lv + 0{e 2 ) 1 
where is the Schwarzschild metric, e is some small parameter and h^ u is the general first order 
perturbation of the metric. We put in the Einstein vacuum equations 5i^„ = 0, where 3?^,, is the 
Ricci tensor, and obtain Si^y = + eR^ + 0(e 2 ) = 0, which leads to the well known Regge- Wheeler 
(RW) differential equations ei?^ t „ = for the perturbation eh^ v . Before trying to solve these equations, it 
is possible to expand e/i M „ in tensor spherical harmonics in the Schwarzschild coordinates (t, R, 9, (f) (see 
Mathews (1962), Zcrilli (1970) and Thorne(1980)). eh^ u falls in one of the two following general classes of 
perturbations, depending on its parity under rotations about the origin performed on the 2-dimensional 
manifold t =const., R =const. : one class (the even-type one) has parity (— 1) and its general (symmetric) 
form is 



ehuv 



I \ 2/R)H { 



im 



\ 



i III 



Hi 

(1 - 2/R)- 1 H* m 
h\ m d B 



h e m de 
h{ m dg 



R 2 (K fm + G lm di) 
R 2 G lm d 2 



h e m d^ 
R 2 G tm d 2 



\ 



im i 



(26) 



R 2 sin 2 9{K lm + G tm d 3 ) ) 



where Ye m (0,<j>) is the standard ^,m-mode spherical harmonic, Hf m (t,R), hf m (t,R), K (t,R) and 
G (t, R) are the corresponding functions of the non-angular coordinates and we define the partial 
differential operators d\ = dgg, d 2 — dg^ — (cos8/sin9)d^, and d 3 = (1/ sin 2 9)8^ + (cos6 / sin9)dg . 

The odd class has parity and its general (symmetric) form is 



where /. 



eh u 



im 



V 






tt m sin6dg 






-f( m (l/sin6)dt 
fi m sin0dg 



(l/2)/| m 5 5 



ti m ain0d e 
f( m sin6dg 
(l/2)/|™3 5 
-f$ m sin 2 9d 4 



Y f , 



(27) 



f[ m {t, R) and 84 = (l/sin9)dg < f ) — (cos9/sin 2 9)d^,, and $5 = (l/sin9)d^,p + cos9dg — sin9dgg. 



The superposition of perturbations are valid in the linearized theory, so in the matrices above we can 
assume the Einstein summation convention in the indices £, m (truncated at the convenience) for the case of 
a more general multi-mode perturbation. We can simplify these matrices further by exploring the freedom 
to make arbitrary, first order in e, coordinate transformations around (t, R, 9, <p). In particular, the matrices 
achieve its most simple or canonical forms in the so called Regge- Wheeler infinitesimal gauge (Regge and 
Wheeler 1957). Nonetheless, we do not need to go to that gauge here. 



If we linearize (n8h in the shell strengths we obtain for the perturbation eh 



/(1-|)P 

(l-%)-\-Q + P) 

R 2 (-Q 

V 



\ 






P) 

R 2 sin 2 9P ) 



(28) 



with P and Q given respectively by (|10|) and (|llj) after dropping the second order terms in the shell 
strengths appearing in those equations (remember that u = R — 1 and v = cosff) . 
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The linear perturbation (|28|) is a compulsory solution of the linearized equations eR^ v — 0, since it 
is the first order term of the expansion of an exact solution of the full equations. Having in mind the 
axial symmetry of our static diagonal perturbation, an inspection of it shows that it is a superposition 
of even-type modes only. Then, it should be fitted by the matrix (p6| ) with the terms £ = 0, 1, 2, 3 and 
m = retained. We find that this fitting is impossible in general. In fact, the linearized solution (p8h is a 
Regge- Wheeler perturbation only if it satisfy the following additional constraints (here k = 1): 

7o = 0, V+l<D = Q. (29) 
5 

On the other hand, the RW formalism lies on a sound mathematical basis, namely the multipolc 
expansion theory (for a review with emphasis in general relativity, see Thorne 1980), so the reason for 
this drawback must have a physical origin. In fact, as it is shown in Sokolov and Starobinskii (1977), 
to conical singularities there correspond a curvature (and hence a Ricci) tensor proportional to a Dirac 
delta term centered at the symmetry axis. From Einstein's equations, this means that there correspond to 
them a certain energy-momentum tensor with this same structure (it does not matter here how exotic they 
should be) and hence we does not have a true vacuum between the core and the shell, as it is supposed 
from the starting in the RW formalism. By comparing ( p9|) with the conditions (|lj) plus ( [TEl) with k = 1 
for Schwarzschild's metric, we see that the former amounts to rule out the conical singularities from the 
linear approximation. In other words, conical singularities survive to the linearization and their presence 
is an obstruction to the application of the RW formalism. So, what is some surprising in all this is that 
conical singularities strutting the model against its own gravity persist after the linearization, contrary to 
the loosely accepted idea that they should lie in the very non-linear realm of general relativity. 



2.4. The Newtonian limit of the relativistic models 

We now consider the Newtonian limit of the model (|l8| ) in the Schwarzschild coordinates (t,R,6,<fr). 
We assume that there exist a region D in the vacuum between the core's horizon and the shell where the 
conditions of weak gravitational field and slow motion of test particles occur. Then, Eintein's equations 
reduce in D to Laplace's equation for the Newtonian potential $, which relates to the metric g M „ only 
through the temporal component as gtt = 1 + (2/c 2 )<i>. The remaining components of the metric are 
irrelevant to this approximation. 

Now, we consistently assume that D is far away from the horizon, where the Schwarzschild coordinates 
approximate to the usual time plus the Euclidean spherical ones (we maintain the notation (t, R, 9, <p) in 
the later approximation). Next, we pass to cylindrical coordinates (t, r, z, (f>) via z — Rcos9, r — RsinO and 



expand the component g tt of (18) to the first order in v , T>, Q and O, to obtain after some manipulation: 



gu = l + (2/c 2 )$ = 1 + 21* - £- ^1 + 2Vz[l - 2 ][i _ i] 

Q(2z 2_ r 2 }[1 _ 2 ][{1 _l_ ) 2__y ] {30) 

.0(2*3-3*r 2 )[l-|][l-i][l-| + ^]. 



The equation above was presented in VL1 (without the dipole) in a rather obscure (yet correct) form 
for the sake of obtaining the Newtonian limit. By assumption, the non-dimensional constants satisfy 
\i/0,V,Q,O\ <C 1; moreover, R satisfies R S> 2 in the region D, so that each square bracket appearing in 
(30) reduces to unity in this approximation. The only place where R itself survives is in the term 2/R of 
(30), just that due to the monopolar core with mass M. The final step is to rewrite the surviving terms 
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with the unit of length L = GM/c 2 appearing explicitly (remember that until now R stood for R/L, etc.). 
The result for <E> is 

$ = c 2 u - ^ + ?fz+ gj(2z 2 - r 2 ) + gj(2z 3 - 3zr 2 ). (31) 

For comparison, let us briefly recall ab initio the proper Newtonian formulation: let the coordinate origin 
stay at the center of mass of the monopolar core (with mass M), z be the symmetry axis of the core-shell 
system and Dm be the region between the smallest and the largest spheres centered at the origin that isolate 
the inner vacuum from the core and the shell. We have to solve Laplace's equation in Dm for the axially 
symmetric Newtonian potential. By using the standard Legendre expansion, we arrive at the following 
gravitational potential $m felt by test particles evolving in Dm- 

C M 1 1 

$Ar = ~ G[/ ° + hZ + 2 h{2z2 ~ ^ + 2 h{2z3 ~ ^ + "' ] > (32) 

where R 2 = r 2 + z 2 = x 2 + y 2 + z 2 , and x,y,z are the usual Cartesian coordinates. Iq, I\, h, I3 are 
respectively the constant, dipole, quadrupole and octopole shell strengths given by the following volume 
integrals over the shell with mass distribution p{R, 0) (P„ standing for the Legendre polynomial of order n): 

I n = JJJ P (R,6)^^ldV. (33) 

shell It 

Obviously, the regions D and Dm above must have a non-empty intersection if the Newtonian approximation 
to the full relativistic case is valid, so we assume this and compare both expressions ( 511 ) and ( p2| ) for the 
potential in D n Dm, thus obtaining 

l^o = -^x J L^|«l, 

|2? = -f xL 2 £|«l, (34) 



L 3 ^\ < 1, 

M ^ 



where m is the mass of the shell. 



The question is: for which Newtonian shells are these constraints on the constants /„ attainable? To 
exemplify, let the shell be an homogeneous ring of mass m, radius a, and centered on the z-axis at z = 6, 
whose density in cylindrical coordinates is p — ^^^(^ — a)d(z — b). The integrals I n in this case are 

Io/m = {a 2 + b 2 )- 1 ' 2 = a- x (l + b 2 /a 2 )~ 1/2 , 
hjm = b(a 2 + b 2 )- 3 / 2 = ba- 3 (l + 6 2 /a 2 )- 3 / 2 , 

him = \{2b 2 - a 2 ){a 2 + b 2 )-^l 2 = ~\a~ 3 {\ - 26 2 /a 2 )(l + 6 2 /a 2 )~ 5 / 2 , ( } 

him = i(26 3 - 3a 2 6)(a 2 + & 2 )~ 7 / 2 = "-\ba^{Z - 26 2 /a 2 )(l + b 2 /a 2 )- 7 ' 2 . 

In the limit a — ► the ring reduces to a particle with mass m placed at z — b, while if b — » it reduces 
to a ring placed at the equatorial plane with vanishing odd multipoles. In the case of m/ M < 1, and 
given the characteristic core lenght L — GM/c 2 (a relativistic parameter we know, but that can be crudely 
anticipated already in the Newtonian frame by assuming light can be trapped by gravity like ordinary 
matter), we see that all constraints in ( j34| ) are satisfied for, say, a^> \b\,L (for the limiting case of a point 
particle at z = b, we let a — + and this condition become |6| S> L). Although intuitively expected, this 
example enlightens and realizes the criteria to overlap both theories, that is to say, to assure the validity 
of the assumption D n Dm ^0, and should be compared to the treatment made for example in Perry and 
Bohun (1992) for Weyl's solutions with usual, decreasing core-type multipoles. 



- 13 - 



We discuss now the physical role played by the constant vq in the model. Obviously, v — z/q, 7 = 70 
in (Q) is a solution of eqs. (Q) and (||) for any values of the constants vq, 70. By the preceding discussion, 
its Newtonian limit leads to the constant potential $n = c 2 v — —GIq. What is the physical meaning of 
a constant Newtonian potential? It does describe either trivially the empty space or rather the interior of 
a very special class of matter distributions, the so called homoeoids (see Chandrasekhar 1987 and Binney 
and Tremaine 1987). The value of the constant potential inside an homoeoid is not arbitrary, having its 
value fixed by continuity requirements of the potential through the mass distribution of the homoeoid. 
Homoeoids are thus gravitationally undetectable from inside. This shows that the constant solution above 
is a struted, relativistic version of a Newtonian homoeoid. If in addition we remove the conical singularity 
from inside the relativistic homoeoid by letting 70 = (remember that conical singularities are not globally 
removable in static relativistic shell solutions since they are indispensable to strut the shell against its own 
gravity, so 70 = really moves them outside the shell), we see that a simple rescaling of the time and the 
radial coordinate reduces the metric inside the homoeoid to that of Minkowski. On the other hand, the 
condition ( pi| ) with T> = Q = O = shows that if we add a black hole inside an homoeoid then the metric 
does not extend analytically through the black hole horizon unless we maintain the "strut" in place with 
strength 70 = 2^o. This is a non-intuitive aspect of the (relativistic) composition of a black hole with the 
rather simple homoeoidal shell. 

An easy point that is worth to emphasize concerning the interpretation of multipolar expansions of 
core and shell types is the following: core multipoles measure deviations from sphericity of central mass 
distributions, while shell multipoles measure how much shells deviate from homoeoids, the later being a 
rather large class of distributions in which homogeneous spherical shells are very particular members. 

We close this section with a brief discussion about the Henon-Heiles structure of the present Newtonian 
limit. It was already pointed out in VL1 that the potential of the shell alone pertains to the Henon-Heiles 
family of potentials, yet it does not suffice from its own to confine orbits in the intermediate vacuum (see 
the opposite signs of r and z in the quadrupole term of ([u|)). Though in practice we need to consider the 
full potential of a given galactic model viz. its effective potential & e ff- The hamiltonian H of a test particle 
with mass [i is given in this case by 

where £ is the conserved angular momentum of the particle associated to the axial symmetry of the galaxy 
(Binney and Tremaine 1987). It is also usual to assume that galaxies have further a reflection symmetry 
about an equatorial plane, in which case the motion restricted to that plane is integrable, in particular 
having a central stable circular orbit at some fixed radius, say, r = rg. To study a nonplanar stelar orbit as a 
small deviation from the planar one, we perform a series expansion of <& e ff in the variables ((r — ro, z), thus 
approximating the full motion by a bi-dimensional harmonic oscillator perturbed by higher order terms. 
This is the very astronomical origin (after truncating the series and idealizing the numeric coefficients) of 
the cubic Henon-Heiles polynomial, now a paradigm of non-integrable potential. Some history about this 
potential in astronomy can be traced, for example, from Contopoulos (1960), Barbanis (1962), van de Hulst 
(1962), Ollongren (1962) and Henon and Heiles (1964). 

By adding the other terms to the terms originated from the shell expansion to form <& e f f , we are able 
to confine test motions in the intermediate vacuum around the central stable orbit. If the galaxy does not 
have reflection symmetry around a middle plane (Ii,l3 7^ in the Newtonian case and T>, O 7^ in the 
relativistic one), planar central stable orbits are not possible, there remaining only stable orbits of distorted, 
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nonplanar type. We will return in section |3| to the discussion of the implications for orbit regularity of the 
largely assumed hypothesis of reflection symmetry about middle plones in galaxy modeling. 



3. Numerical results: reflection symmetry and chaotic motion in relativistic and Newtonian 

core-shell models 

We study the chaoticity of orbits in the intermediate vacuum of Newtonian static axially symmetric 
core-shell models and then compare them to the geodesies of the corresponding relativistic cases. It is 
amazing that this class of models has relativistic and Newtonian counterparts to be compared, there existing 
moreover sound observational motivations for both. 

Specifically, we analyse in both cases the role played by a discrete symmetry, namely the reflection 
symmetry around the equatorial plane, in regularizing orbits against chaoticity. The hypothesis of reflection 
symmetry is a widely spread assumption in astronomical modeling and its twofold justification lies firstly 
in that too many real objects to be modeled seem in fact nearly symmetric with respect to a plane. On 
the other hand, the models are simplified in both analytical and numerical aspects by that assumption. 
Of course, we do not expect that symmetry be realized in nature exactly so it is relevant to search for 
possible orbital effects of its breaking. It is important to realize that all cases treated here preserve the 
axial symmetry — the next continuous symmetry after the missing spherical one, which allows us to isolate 
the dynamical effects of only breaking or preserving the reflection symmetry itself. 

Another advantage of the present study is that it is model independent (of course within the class of 
models we are dealing with) as we are focusing on general multipolar shell expansions instead on specific 
mass distributions. 

All is made non-dimensional again in this section by formally taking $n <-> $jv/c 2 ,i? «-> R/L, etc., 
H <-> H/(iic 2 ),p z «-> p z / (pc), etc., and i «-+ £/ (pcL) in equations ([3l]-32) plus (|36|), such that the Newtonian 
system we consider comes from the hamiltonian (R 2 = r 2 + z 2 ) 



H=i(p 2 r +p 2 z ) + l^-l 

+ Vz + ±Q(2z 2 - r 2 ) + \0(2z z - 3zr 2 ). 

On the other hand, the geodesic system comes from the lagrangean £ <-> C/ (pc) given by 



(37) 



£ = \g a pi a i^ (38) 

where the metric tensor g a p and the coordinates x° — t^x 1 = p,x 2 — z,x 3 — tfi are obtained from the 
non-dimensional invariant interval ([!]), and the dot stands for the derivative d/ds. The proper definition 
of time- like geodesies out of the metric interval furnishes the first constant of motion £ = 1/2. The 
Euler-Lagrange formulation of the geodesic equations is 

d0C £-0, (39) 



ds dx^ dx^ 

from which the two additional constants of motion associated to the static nature and axiality of the 
relativistic system are read out, namely the relativistic energy h and angular momentum I defined by 



(40) 
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The remaining two equations in ( |39| ) describe the dynamics for the variables p, z. We note that these p, z 
are Weyl's coordinates. From (^) we see that whenever the Schwarzschild coordinate R satisfies R ^> 2 
it approximates the usual radial spherical one (denoted here by the same letter R) and the Weyl p, z 
approximate the usual cylindrical ones r, z. We remember also that in the Newtonian limit dt/ds sa l/^/g^ 
(low velocities), gu ~ 1 + 2$n/c 2 , so the relativistic energy h and the Newtonian energy E «-» /i$Ar/(/zc 2 ) 
of an orbiting particle with mass p are related through h 1 + E, and the corresponding angular momenta 
through I ss ^. We put = and k = 1 in this section. We impose the conditions ( |l4| ) and ([l5|) for the 
absence of conical singularities on all running relativistic situations presented bellow, except of course for 
the purely dipolar shell. All Newtonian as well as relativistic Poincare sections (surfaces of section) shown 
here are made at the plane z = (with the appropriate coordinate interpretation in each case). 

In Fig. 1 we present typical effective Newtonian potential wells for the bounded orbits we are considering 
{U stands for $ e // and the positive parts of the potential surfaces surrounding the wells are cut). Figs. 1(a) 
and (c) show cases of oblate and prolate potentials, respectively, both possessing reflection symmetry about 
the equatorial plane given by z — 0. Oblate cases (a) have positive quadrupole strengths (Q > 0) while 
prolate cases (c) have Q < (respectively 1% < and I2 > 0, see eqs. ( |3l| ) and ( p2[ ) and, for the specific 
case of a ring, eq. (]3^)). The presence of the reflection symmetry implies the vanishing of all odd shell 
multipoles (T>, O and I\, I3 here). We see that the oblate case has only one unstable equilibrium point (in 
the equatorial plane) while the prolate case has two of them (symmetrically placed outside the equatorial 
plane), which is easy to understand in terms of the exterior oblate/prolate shell mass distributions. Figs. 
1 (b) and (d) show the typical deformation of the previous oblate and prolate cases, respectively, when we 
break the symmetry of reflection about z = by introducing non- vanishing odd shell multipoles (2?, O). 

On the other hand, the relativistic regime is felt in the situation we are concerned with mainly through 
the existence of one more unstable equatorial equilibrium point in addition to those already present in the 
Newtonian potential. This intrinsically relativistic additional unstable point is associated to the presence 
of the black hole at the center and marks the point above which the orbits fall into the black hole. In the 
remaining, the relativistic potentials are qualitatively similar to those displayed in Fig. 1 for the Newtonian 
case. 

We have made a wider exploration of the shell parameters, energy and angular momentum than shown 
here, with the same conclusions in large. The parameter values we actually chose to show are a compromise 
between the need to make consistent comparisons and the sharpening of the effects we found. 

The case of a purely dipolar shell has been anticipated in VL2, where we stressed that the monopolar 
core plus shell dipole is non-trivial already in the Newtonian context. Moreover, we show that the 
Newtonian case is integrable whereas the relativistic one is chaotic, which justifies the characterization of 
chaos in the purely dipolar case as an intrinsic general relativistic effect. Here, we explore the fact that 
the Newtonian dipolar shell breaks the reflection symmetry of the model without breaking the integrability 
of the motion itself. Thus core-shell models provide us with two very distinct integrable situations from 
the point of view of the reflection symmetry: core plus dipolar shells (D-cases), which do not have that 
symmetry, and purely monopolar cores (Kleperian cases) which do have it. A typical Poincare section of an 
integrable D-case is shown in Fig. 2(a). 

Firstly, we perturb both integrable configurations above with a reflection symmetry preserving, oblate 
quadrupolar term (Q > 0). Surfaces of section for core plus oblate quadrupolar shells (oblate Q-cases) are 
shown in Figs. 2(b,c,f) while sections for core plus oblate dipolar-quadrupolar shells (oblate DQ-cases) are 
shown in Figs. 2(d,e). These figures show that the breaking of the integrability of the D-cases (without 



- 16 - 



reflection symmetry) by a quadrupole is much stronger than that of the Keplerian case (with reflection 
symmetry) by the same quadrupole. Moreover: Fig. 2(c) shows that chaoticity in the oblate Q-cases is 
in fact present only in a residual, "microscopic" level, Fig. 2(f) illustrates the robustness of the strong 
regularity of the oblate Q-cases against the varying of the quadrupole strength, while Fig. 2(e) shows that 
the strong chaoticity in the DQ-cases is very dependent on the energy, in contrast with the robustness of 
the almost regularity of the oblate Q-cases also with the energy. 

Fig. 3 shows the effects of introducing a reflection symmetry breaking, octopolar perturbation on the 
two integrable configurations (D-cases and Keplerian cases), and also on the almost regular oblate Q-cases. 
Figs. 3(a) and (f) show the effects of a purely octopolar shell on the Keplerian case (O-case) at two different 
octopole strengths. We see that these effects are very much stronger in comparison with the almost-regular 
oblate Q-cases. Figs. 3(b-e) show the effects of the various combinations of shell multipoles and confirm 
the dominance of the chaotic effects associated to the odd multipoles over those related to the purely 
oblate Q-cases. In particular, the enhancement of chaoticity associated to the breakdown of the reflection 
symmetry is reinforced by the oblate QO-cases presented in Fig. 3(c) in comparison with the oblate Q-cases, 
and Figs. 3(d) and (e) confirm the strong dependence with the energy of the chaoticity associated to the 
lack of reflection symmetry. 

Fig. 4 shows surfaces of section for prolate quadrupolar cases (Q < 0). Figs. 4(a-c) show the very much 
higher chaoticity of the prolate Q-cases in comparison to the previous quasi-regular oblate Q-cases. This is 
to be expected in view of the prolate Q-cases having two unstable equilibrium points (and hence suffering 
the simultaneous influence of the two corresponding instability regions) instead of only one occurring in the 
oblate Q-cases. These figures also show that the central primary stable orbit, deeply inside the accessible 
region, is strikingly the first one to bifurcate and drive the chaotic behavior with the rising of the energy. 
The strong dependence with the energy of the whole set of bifurcations is also a distinctive aspect of 
these figures. In Fig. 4(d) we show the effect of applying the perturbing prolate quadrupolar term on the 
integrable, reflection symmetry missing D-cases (prolate DQ-cases). We find that the absence of mirror 
symmetry in the prolate DQ-cases causes the opposite to that occurring in the oblate DQ-cases, namely 
the missing mirror symmetry does greatly enhance the orbit regularity in the prolate cases. Fig. 4(e) 
does confirm that the breakdown of mirror symmetry associated to the addition of an octopolar term to 
the prolate Q-cases (prolate QO-cases) strongly reduces the chaotic manifestation of the prolate Q-cases 
themselves. In particular, Figs. 4(d-f) show evidences that the lack of symmetry introduces a robust family 
of regular asymmetric orbits around a stable primary one in the bounded region. This regularizing effect 
due to the lack of mirror symmetry in the prolate DQO-cases is easily understandable in terms of the 
progressive lack of influence of one among the two unstable equilibrium points of the potential on the 
bounded region accessible to the particle if odd shell multipoles are put in scene (cf. Figs. 1(c) and (d)). 

We present in Fig. 5 the findings for the corresponding relativistic core-shell models. We choose for 
these figures the same non-dimensional multipole shell strengths and angular momenta of the corresponding 
Newtonian figures and, except for the purely dipolar shell, we always eliminate the conical singularities (or 
"struts") from the intermediate vacuum by fixing 70 and T> in terms of Q and O through the conditions 
( |l4| ) and (|lj). Also, we always put vq = here. The relativistic energy h is chosen in each case to be close 
that given by the approximate relation h w 1 + E where E is the Newtonian energy. 

Fig. 5(a) shows a typical surface of section for a black hole plus purely dipolar shells (RD-cases), 
its strong chaoticity having to be compared with both the almost regularity of the oblate relativistic 
quadrupolar case (RQ-case) shown in Fig. 5(b) (its very small chaotic zones are visible if we zoom the 
cross-type regions between the islands) and the corresponding integrable Newtonian D-case shown in Fig. 
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2(a). [In particular, the smallness of the chaotic zones of the oblate RQ-case shown in VL1 occasioned the 
misleading statement made by us therein of these cases not exhibiting chaos, which was promptly corrected 
in Vieira and Letelier 1996b.] 

In fact, Fig. 5(b) does not represent the full chaotic behavior possible to oblate RQ-cases. As we saw 
above about the relativistic counterparts of the Newtonian potentials, the relativistic oblate cases have one 
more unstable scape point (in fact, an infall point into the black hole) inermost at the equator in addition 
to that one provided in both cases by the oblate shell itself. Hence, if we vary the energy and the angular 
momentum so that both unstable regions are equally accessible by the bound nearly equatorial geodesies, it 
is possible to enhance strongly the chaotic behavior of Fig. 5(b) to encompass even a significant portion of 
the more external region of the surface of section (that associated to the nearly equatorial orbits). We will 
see more about this enhancement for the relativistic cases below in Fig. 5(e). 

Fig. 5(c) shows a surface of section for relativistic dipolar-octopolar cases (RDO-cases). Remember 
that both odd components are simultaneously needed in view of the off-strut constraint (|l^). Fig. 5(d) 
is the same for full oblate RDQO-cases. Figs. 5(a-d) show that the chaotic zones are greatly enlarged 
whenever odd shell multipoles are present, in contrast with the almost-regular oblate RQ-case shown. 

Figs. 5(e) shows a surface of section for relativistic prolate cases (prolate RQ-cases). It is to be 
compared to the corresponding prolate Newtonian one (see Fig. 4(b)), in particular we note the evident 
fingerprints of the bifurcation series starting from the originally stable central primary orbit already present 
in the prolate Q-cases. We also see that the prolate RQ-cases share with the prolate Q-cases the same 
relatively wild chaotic behavior as compared to the respective oblate (R)Q-cases. 

Like in the Fig. 5(b), the dynamics would be more chaotic in Fig. 5(e) if the angular momentum and 
the energy (unaltered here for the sake of full comparison) were ajusted to allow all the three unstable 
points of the relativistic potential (instead only the two ones related to the prolate character of the shell, see 
Fig. 1 and the comments therein about the additional black hole unstable point) have stronger combined 
influence on the motion. So, the remaining differences between the relativistic figure 5(e) and the Newtonian 
one 4(b) are due to the following : the central region (in the surface of section) is the only one able to 
become chaotic in the Newtonian case since its orbits are just those (non-planar ones) that can reach the 
unstable regions near the two unstable scape points symmetrically placed outside the equatorial plane, 
while the more external (in the surface of section) quasi-equatorial orbits are regular since they are far 
from those two points and under the strong influence of the equatorial stable orbit in the boundary of the 
surface of section. When we change the Newtonian central mass by a true black hole then a third unstable 
scape point (in fact a rolling down point into the black hole) is placed innermost in the equatorial plane. 
Associated to it there is one more unstable region which can be reached, this time by the quasi-equatorial 
orbits. The specific combination of parameters of Fig. 5(e) is such that the central, non-planar orbits can 
feel only moderately the influence of the two scape points linked to the presence of the prolate shell, while 
the quasi-equatorial orbits are feeling strongly the nearby influence of the unstable point associated to the 
black hole. As we said above, other combinations of parameters will allow a much stronger spread and 
eventually the overlaping of both types of chaotic regions. 

Finally, Fig. 5(f) shows the restoration of a robust set of regular asymmetric orbits around a primary 
stable one when we break the reflection symmetry around the equatorial plane by switching on the odd 
multipoles onto the prolate RQ-cases, much the same as occured in the prolate Newtonian Q-cases. 
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4. Discussion 

In the first part of this article we presented an unifying discussion concerning with the properties 
and, more important, the physical content of some relativistic, static, axially symmetric core-shell models 
on their own as well as in connection with the corresponding linearized and Newtonian models taken as 
limiting cases. The analytical and observational motivations for these relativistic and Newtonian models 
were also shown. This was acomplished in a reasonably self-contained manner in the introduction and 
section 2 and needs no additional comments. 

In the second, numerical part of this work we explored the fact that the models i) are generic within 
the class they pertain and ii) have exact relativistic and Newtonian counterparts to be compared, to study 
the chaoticity of bound orbits in the vacuum between the core and the shell. 

Spccificaly, we firstly tested the relevance of the presence/absence of the reflection symmetry around 
the equatorial plane for the chaotic behavior of the orbits. We find consistent evidences for a non-trivial 
role played by the reflection symmetry on the chaoticity of the dynamics in both the relativistic and 
Newtonian cases. We summarize these findings as follows: the breakdown of the reflection symmetry about 
the equatorial plane in both Newtonian and relativistic core-shell models docs i) enhance in a significant 
way the chaoticity of orbits in reflection symmetric oblate shell models and ii) inhibit significantly also the 
occurrence of chaos in reflection symmetric prolate shell models. In particular, the lack of the reflection 
symmetry provides the phase space in the prolate case with a robust family of regular orbits around a 
stable periodic orbit that is otherwise missing at higher energies. 

The other point we addressed about the chaotic behavior of orbits was for the consequences of 
substituting true central black holes by Newtonian central masses. We find that the relative extents of 
the chaotic regions in the relativistic cases are significantly larger than in the corresponding Newtonian 
ones. Although not surprising in thesis, the strong differences between both regimes are in order in view 
of the procedure found in the literature of simulating the presence of a black hole at the core of a galaxy 
with a naive — 1/r Newtonian term (see for example Gerhard and Binney 1985, Sridhar and Touma 1997 
and Valluri and Merrit 1998). This approximation is certainly valid far from the black hole but not in its 
proximity as is the case here. 

These findings stress i) the non-trivial role of the reflection simmetry in both relativistic and Newtonian 
regimes, in contrast with its universal acceptance in astronomical modeling, ii) the strong qualitative 
and quantitative differences between relativistic and Newtonian regimes, in particular when dealing with 
orbits in the vicinity of black holes, and iii) the intrincate interplay between both aspects when they are 
simultaneously present. 

The true dynamical aspect related to the role of reflection symmetries on the regularity of the orbits 
refers of course to the parity of the constants of motion (the famous third one in the present axially 
symmetric case) with respect to the coordinates, as well as their power to regularize orbits in phase space. 
Translated to these terms, the findings above are saying that the additional constant of motion in question 
is by far more powerful whenever i) it is an even function of the coordinates for oblate cases, or ii) it is 
an odd function of the coordinates for prolate cases. Since the issue of finding global or even approximate 
constants of motion given the dynamics is hard if not feasible in most cases, the former procedure of simply 
checking the presence/absence of reflection symmetries is of much more practical interest. 

Additional study is needed to see if and how our findings are extendible to more realistic configurations, 
with and without black holes/central masses. The obvious improvement is to fill the intermediate vacuum 
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with some reasonable mass distribution. While this is far from obvious in general relativity, it is easy to 
do in the Newtonian context. For example, we are considering just superpose our shell multipoles to some 
relevant potentials of celestial mechanics such as Plummcr-Kuzmin, Ferrers and others and even triaxial 
potentials. They have one or more planes of symmetry, with or without axial symmetry, and so we can 
perturb them with shell terms to verify similar effects to those we found. 

Here we considered only tube orbits (l z ^ 0) since we are interested in the effects of the existence 
of planes of symmetry on tridimensional motion (somewhat similar effects to the found here are seen in 
Gerhard 1985 for the restricted case of planar motions with respect to the existence of one or more lines 
of symmetry) . The model improvements above will allow us to study the effects of breaking the reflection 
symmetry also on tridimensional box orbits. 

Meantime, the fully relativistic program is in progress. We have recently succeeded in giving rotation 
to a black hole (i. e. in converting it into a Kerr black hole) plus a dipolar shell term (Letelier and Vieira 
1997). There are increasing evidences for the existence of black holes, particularly inside active galactic 
nuclei (Kormendy and Richstone 1995), which motivate us to consider also rotating core-shell models in the 
same lines followed here. 

The authors thank FAPESP and CNPq for financial support, and Jorge E. Horvath and Andre L. B. 
Ribeiro for discussions. 



REFERENCES 

Araujo, M. E., Letelier, P. S., and Oliveira, S. R., 1998, Class. Quantum Grav., in press. 
Arnaboldi, M., Capaccioli, M., Cappellaro, E., Held, E. V., and Sparke, L., 1993, A&A, 267, 21. 
Arnol'd, V. I., 1963a, Russ. Math. Surv., 18, 9. 
Arnol'd, V. I., 1963b, Russ. Math. Surv., 18, 85. 
Barbanis, B., 1962, Zeit. Astrophys., 56, 56. 

Barlow, M. J., Drew, J. E., Meaburn, J., and Massey, R. M., 1994, MNRAS, 268, L29. 
Barnes, J. E., 1993, ApJ, 419, L17. 

Barnes, J. E., and Hcrnquist, L., 1992, ARA&A, 30, 705. 

Berry, M. V., 1978, Am. Inst. Phys. Conf. Proc, 46, 16. 

Binncy, J., 1982a, MNRAS, 201, 15. 

Binncy, J., 1982b, MNRAS, 201, 1. 

Binncy, J., and Spergcl, D., 1982, ApJ, 252, 308. 

Binney, J., and Tremaine, S., 1987, Galactic Dynamics (Princeton: Princeton Univ. Press). 
Carminati, J., and McLcnaghan, R. G., 1991, J. Math. Phys., 32, 3135. 
Chandrasekhar, S., 1987, Ellipsoidal Figures of Equilibrium (New York: Dover). 
Chevalier, R. A., 1997, Science, 276, 1374. 
Contopoulos, G., 1960, Zeit. Astrophys., 49, 273. 
de Zeeuw, T., 1985, MNRAS, 216, 273. 



-20- 



de Zeeuw, T., Peletier, R., and Franx, M., 1986, MNRAS, 221, 1001. 

Dupraz, C, and Combes, F., 1987, A&A, 185, LI. 

Evans, N. W., Halner, R. M., and de Zeeuw, 1997, MNRAS, 286, 315. 

Gerhard, O. E., 1985, A&A, 151, 279. 

Gerhard, O. E., and Binney, J., 1985, MNRAS, 216, 467. 

Groenewegen, M. A. T., Whitclock, P. A., Smith, C. H., and Kerschbaum, F., 1998, MNRAS, 293, 18. 

Henon, M., and Hciles, C, 1964, AJ, 69, 73. 

Kolmogorov, A. N., 1954, Dokl. Akad. Nauk. SSSR, 98, 527. 

Kormendy, J., and Richstonc, D., 1995, ARA&A, 33, 581. 

Madsen, J., 1991, ApJ, 367, 507. 

Moon, P., and Spencer, D. E., 1988, Field Theory Handbook (Berlin: Springer- Verlag) . 

Letcher, P. S., and Vieira, W. M., 1997, Phys. Rev. D, 56, 8095. 

Lctelier, P. S., and Olivcira, S. R., 1998, Class. Quantum Grav., 15, 421. 

Lyndcn-Bcll, 1962, MNRAS, 124, 95. 

Malin, D. F., and Carter, D., 1983, ApJ, 274, 534. 

Mathews, J., 1962, J. Soc. Indust. Appl. Math., 10, 768. 

Merrit, D., and Fridman, T., 1996, ApJ, 460, 136. 

Mcrrit, D., 1997, ApJ, 486, 102. 

Meyer, F., 1997, MNRAS, 285, Lll. 

Moser, J., 1967, Math. Ann., 169, 136. 

Norman, C. A., Sellwood, J. A., and Hasan, Hashima, 1996, ApJ, 462, 114. 
Ollongrcn, A., 1962, Bull. Astron. Inst. Netherlands, 16, 241. 

Panagia, N., Scuderi, S., Gilmozzi, R., Challis, P. M., Garnavich, P. M., and Kirshncr, R. P., 1996, ApJ, 
459, L17. 

Perry, G. P., and Bohun, C. S., 1992, Phys. Rev. D, 46, 1866. 

Poincare, H., 1957, Les Methodes Nouvelles de la Mecaniquc Celeste (New York: Dover), vols. 1-3. 

Quinn, P. J., 1984, ApJ, 279, 596. 

Ralston, J. P., and Smith, L. L., 1991, ApJ, 367, 54. 

Regge, T., and Wheeler, J. A., 1957, Phys. Rev., 108, 1063. 

Reina, C, and Treves, A., 1976, Gen. Rel. Grav., 7, 817. 

Reshetnikov, V., and Sotnikova, N., 1997, A&A, 325, 933. 

Richstone, D. O., 1982, ApJ, 252, 496. 

Robertson, H., and Noonan, T., 1968, Relativity and Cosmology (London: Saunders). 
Sackett, P. D., and Sparke, L. S., 1990, ApJ, 361, 408. 
Schwarzschild, M., 1979, ApJ, 232, 236. 
Schwarzschild, M., 1993, ApJ, 409, 563. 



-21 - 



Sokolov, D. D., and Starobinskii, A. A., 1977, Sov. Phys. Dokl., 22, 312. 

Sridhar, S., and Touma, J., 1997, MNRAS, 287, LI. 

Suen, Wai-Mo, 1986a, Phys. Rev. D, 34, 3617. 

Suen, Wai-Mo, 1986b, Phys. Rev. D, 34, 3633. 

Thome, K. S., 1980, Rev. Mod. Phys., 52, 299. 

Valluri, M., and Merrit, D., 1998, ApJ, to appear. 

van de Hulst, H. C, 1962, Bull. Astron. Inst. Netherlands, 16, 235. 

Vieira, W. M., and Letelier, P. S., 1996a, Phys. Rev. Lett., 76, 1409 (VL1). 

Vieira, W. M., and Letelier, P. S., 1996b, Los Alamos Data Bank flgr-qc/9608030| ). 

Vieira, W. M., and Letelier, P. S., 1997, Phys. Lett. A, 228, 22 (VL2). 

Vishveshwara, C. V., 1970, Phys. Rev. D, 1, 2870. 

Zcrilli, F. J., 1970, Phys. Rev. D, 2, 2141. 

Zhang, X.-H., 1986, Phys. Rev. D, 34, 991. 



This preprint was prepared with the AAS IATjrjX macros v4.0. 



-22 - 



Fig. 1. — Typical effective Newtonian potential wells we are dealing with. U stands for $ e ff and the positive 
parts of the full potential surfaces are cut. (a) shows an even (T> = O = 0) oblate (Q > 0) shell and (b) 
shows its perturbation due to the presence of the odd shell multipoles T> and/or O. On the other hand, (c) 
shows an even prolate (Q < 0) shell and (d) shows its perturbation due to the presence of the odd shell 
multipoles. Note that (a) and (c) have a plane of symmetry, which is broken by the odd multipoles in (b) 
and (d). 

Fig. 2. — This figure exhibits, for the parameter values shown, surfaces of section at the plane z = of 
integrable Newtonian configurations (D-cases and Keplerian cases) perturbed by oblate shell quadrupoles 
only. Note that quadrupolar perturbations do preserve the mirror symmetry. In all figures L accounts for 
the non-dimensional angular momentum I. (a) shows a typical section of the integrable D-case. This case 
does not have a plane of symmetry, in contrast with the also integrable spherically symmetric Keplerian case, 
(b), (c) and (f) show the perturbed Keplerian cases (oblate Q-cases), while (d,e) show perturbed D-cases 
(DQ-cases). Note that orbit regularity is strongly broken (preserved) in the absence (presence) of mirror 
symmetry, this conclusion being robust against multipole strength variations. 

Fig. 3. — In this figure we add octopolar components (O ^ 0) to the shell, which necessarily does break 
the mirror symmetry. In (a) and (f) we see the strong chaotic effect of perturbing Keplerian cases with 
octopoles (O-cases) (compare with the ahnost-regularity of Figs. l(b,c,f)). In (b) we see the strong chaotic 
effect of octopoles on integrable D-cases (DO-cases), which is fully expected as mirror symmetry is already 
absent from the starting. In (c) we break the reflection symmetry of the almost-regular oblate Q-case with 
octopoles (QO-cases) and in (d,e) we present the full DQO-case for two slightly different energies, which 
confirms the strong dependence of the chaoticity with the energy when the mirror symmetry is broken, in 
contrast with the robustness of the almost-regular, mirror symmetric cases against energy variations. 

Fig. 4. — We show in this figure prolate Q-cases. In (a-c) we illustrates the well known fact that (mirror 
symmetric) prolate Q-cases are strongly chaotic on their own and that their chaoticity is highly energy 
dependent. We note in particular that the orbit bifurcations toward the chaotic behavior unusually start 
from the central primary stable orbit itself rather than from the boundary. In (d-f) we break the mirror 
symmetry of the prolate Q-cases with different combinations of odd multipoles. This does cause a strong 
regularizing effect on the orbits, in particular by restoring, of course in a distorted fashion, the whole family of 
regular orbits around a primary stable orbit. Moreover, this restoration is robust against multipole strength 
as well as energy variations. 

Fig. 5. — We present some surfaces of section at the plane z = for the corresponding fully relativistic 
core-shell configurations. In this figure, the coordinates r, z stand for the Weyl ones p, z (see text for more 
details). In addition, in (a) we set vq = 70 = (in the purely dipolar case conical singularities (CSs) are 
unavoidable), while in (b-f) we set v$ = and, in order to eliminate CSs, 70 and T> are given in terms of Q and 
O according to the conditions ( [Til) and (|l5|). The relativistic case confirm the role of the presence/absence of 
the reflection symmetry on the chaoticity of the orbits already detected in the preceding Newtonian figures. 



